ADA014558 


NPS55SY75071  " / 

NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


» 


I 


N 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  Isham  Linder  Jack  R.  Borstlng 

Superintendent  Provost 


Reproduction  of  all  or  part  of  this  report  Is  authorized. 
This  report  was  prepared  by: 


Bruno  0.  Shubert,  Associate  Professor 
Department  of  Operations  Research  and 
Administrative  Sciences 


Reviewed  by: 




rid  A.  Schrady,  Chalrmai 


Department  of  Operations 
and  Administrative  Scl 


search 

ces 


Released  by: 


Robert  Possum 
Dean  of  Research 


UNCLASSIFIED 

security  classification  of  this  page  rlFhwi  Pm«  Enuraip 


1-^ 


REPORT  DOCUMENTATION  PAGE 


iKEJA 


] NP^5SY75071 


[2.  GOVT  ACCESSION  NO 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1 RECIPIENT'S  CAT ALOG  NUMBER 


I "TiTLe  ButNinr 


MODELING  A RANDOM  SEARCH,  / 


S.  TJCilE  BP  fH?PBWT  S PBWOB  CpVEREO 

'echnlcal 


’ NUMBER 


7.  AUTMORriJ 


B.  CONTRACT  OR  GRANT  NUMBERr*; 


a 


Bruno  O.y^hubert  j 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Dept,  of  OR/AS  Code  55Sy 

Naval  Postgraduate  School 
Monterey,  CA  93940 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  S WORK  UNIT  NUMBERS 


£ 


I I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


rbWt: 


Ju 


75 


R OF  PAdES 


25 


14.  monitoring  agency  name  « AOORESSrU  <«/f*rMir  from  Conlrelllnf  Olllem) 


IS.  security  class.  (oI  Ihli  ropertj 

Unclassified 


IS*.  DECLASSI  FI  CATION/ downgrading 
SCHEDULE 


l«.  DISTRIBUTION  STATEMENT  (ol  Ihli  froporO 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  Iho  tbHrmet  miond  In  Block  30,  II  dlllmroni  from  Roport) 


18.  SUPPLEMENTARY  NOTES 


'71 


9.  KEY  WORDS  (Conllnuo  oi 

Random  Search 
Wiener  Process 


•fd«  U n«e«««ary  md  Id9ntlty  by  btoek  numbmr) 


20l  ABSTPIACT  (Continue  on  rovorvo  //  and  Idonfl^  by  6locA  numbmr) 

^A  model  for  random  search  Is  considered  which  differs  from  the 
classical  random  search  by  modeling  the  random  paths  of  the  searcher 
and/or  target  as  those  of  a Wtfiner  process.  The  quantity  of  Interest 
Is  the  distribution  of  detection  time  for  a cookie  cutter  mode  of 
detection.  Only  one-dlmenslonal  search  Is  considered  here.,^ 


DD  ,:?h7,  1473  EDITION  OF  I NOV  68  IS  OBSOLETE 

S/N  0102-014-6601  I 


UNCLASSIFIED 


4.  /'/ 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  CBTiwi  Dmim  Knldtod) 


V. 


1.  Introduction 


This  report  is  a preliminary  study  of  some  models  for  random 
search.  It  begins  by  reexamining  the,  by  now,  classical  formula  of 
Koopman  for  the  probability  of  detection.  It  is  shown  that  the  assump- 
tions under  which  the  formula  was  derived  are  incompatible  with  the 
natural  requirement  that  the  searcher’s  trajectory  be  continuous.  To 
remedy  the  situation  a different  model  of  random  motion  is  needed.  The 
main  theme  of  this  report  is  to  model  such  a motion  as  a continuous 
version  of  a random  walk  - the  Wiener  process.  With  this  modification 
a new  formula  for  the  detection  time  distribution  is  obtained.  This  is 
done  both  for  a stationary  target  and,  under  some  simplifying  assumptions, 
also  for  a moving  target. 

This  work  is  preliminary  in  the  sense  that  only  one-dimensional 
search  is  considered.  Although  search  in  a narrow  straight  may  be  re- 
garded as  one-dimensional,  this  restriction  clearly  limits  the  application 
of  our  model  at  its  present  state.  The  model  Itself,  however,  readily 
extends  to  two  dimensions;  it  is  its  mathematical  analysis  that  becomes 
more  involved.  Thus,  such  a restriction  seemed  necessary  for  the  initial 


phase  of  research. 


r 
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2 . Scatlonary  Target 

In  this  section  we  Intend  to  reconsider  the  basic  formula  for 
random  search 

Wjj 

P ■ 1 - e ^ , (1) 


where  p Is  the  probability  of  detection,  A Is  the  area  in  which  the 
target  Is  located,  W is  the  effective  search  width  and  L is  the 
total  length  of  the  observer's  path.  This  formula  was  originally  derived 
by  Koopman  under  the  following  assumptions  [1]  , page  28: 

(Al)  "The  target's  position  is  uniformly  distributed  in 
A." 

(A2)  "The  observer's  path  is  random  in  A in  the  sense 
that  it  can  be  thought  of  as  having  its  different 
(not  too  near)  portions  placed  independently  of  one 
another  in  A." 

(A3)  "On  any  portion  of  the  path  which  1s  small  rela- 
tively to  the  total  length  of  path  but  decidedly 
larger  than  the  range  of  possible  detection,  the 
observer  always  detects  the  target  within  the 
lateral  range  W/2  on  either  side  of  the  path 
and  never  beyond . " 

Koopman  then  proceeds  by  first  dividing  the  observer's  path  into 
n segments  of  length  L/n  and  using  (Al)  and  (A3)  to  conclude  that  the 
probability  of  detecting  the  target  in  a particular  segment  is  LW/nA. 

He  then  employs  (A2)  to  reason  that  for  sufficiently  large  n the  events 
{detection  in  k-th  segment},  k = 1,  2,  ...  should  be  independent, 
whence  the  formula  immediately  follows  by  letting  n increase  to  infinity. 


WL 


n-x» 


Let  us  now  reexamine  the  assumptions  In  a little  more  detail. 
Assumption  (Al)  seems  quite  reasonable  Inasmuch  as  It  expresses 
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a complete  lack  of  prior  knowledge  of  the  target  location  (Principle  of 
Insufficient  Reason).  However,  It  Is  not  clear  wheth*>r  the  target  is 
stationary  during  the  search  or  whether  It  moves.  In  the  latter  case, 
what  should  be  assumed  about  Its  motion?  Ue  will  return  to  this  question 
shortly. 

Assumption  (A3)  concerns  a mode  of  detection  - the  so-called 
cookie  cutter  detection  - and  can  be  rephrased  more  simply  as: 

(A3')  Detection  occurs  as  soon  as  the  distance  between 
the  target  and  the  observer  decreases  to  or  below 
some  positive,  constant  C = W/2. 

Clearly,  this  is  a matter  of  postulating  a mode  of  detection,  (A3)  or 

'A3')  representing  the  simplest  case. 

Assumption  (A2)  Is  the  crucial  one.  Inasmuch  as  it  defines  the 

random  path  of  the  observer.  Unfortunately,  it  is  stated  rather  vaguely. 

To  try  to  understand  it,  let  us  denote  by  Y^,  t ^ 0 the  observer's 

location  at  time  t,  the  time  t = 0 being  the  beginning  of  search. 

(A2)  seems  to  say  that  and  should  be  independent,  at  least 

for  larger  t > 0.  If  the  observer  moves  with  a constant  speed  v, 

dividing  its  path  into  n segments  of  length  L/n  is  then  equivalent 

to  dividing  the  time  interval  [0,T=L/v]  into  n equal  subintervals  of 

length  T,  n.  Calling  tj^  u point  in  the  k-th  subinterval,  let  us  say 

T 1 

the  midpoint  t,  * — (k+^)  , assumption  (A2)  requires  that  the  random 

variables  Y , k = 0,  . . . , n - 1 be  independent  as  long  as  T/n  is 
‘■k 

not  too  small  (cf.  "not  too  near  portions"  in  (A2)). 


A 


Let  X^,  t 0 be  the  location  of  the  target  at  time  t.  Now, 
if  either 


X.  - X-  ~ uniform  in  A (stationary  target) 


or 


X ,k.-0,  ...,n-l,  independent,  identically  distributed  uniformly 
in  A (moving  target) , 

then  using  (A3) 


P(no  detection  along  a path  of  length  L) 
= P(  I 1 Yj.  - Xj.  I I > I for  all  0 ^ t i T) 
~ P(X  is  for  k - 0 n - 1) 

tk  K 


rectangle  of  length  L/n,  width  W and  center  at  . Notice  that 

*^k 

the  approximate  equality  above  holds  by  (A3)  for  large  n,  i.e.,  for 

small  T/n.  On  the  other  hand,  the  Y 's  are  assumed  Independent 

*"k 

provided  T/n  is  not  too  small.  But  this  prevents  us  from  taking  the 
limit  n -►  “>!  Thus,  we  are  forced  to  drop  the  restriction  on  "not  too 
near  portions"  from  (A2)  and  assume  that  Y^.  and  are  independent 

for  all  T > 0,  no  matter  how  small. 

But,  as  is  well  known,  this  is  incompatible  with  the  obvious 
requirement  that  the  observer’s  path  be  continuous.  Hence,  if  we  indeed 
wish  to  model  a search  for  a stationary  or  moving  target  by  a moving 
craft,  we  have  no  choice  but  to  abandon  (A2)  altogether  and  define  the 
observer's  random  motion  in  a different  way. 
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Again,  let  Y^,  t ^ 0 be  the  location  of  the  observer  at  time  t. 
Since  the  motion  is  to  be  random,  {Y^.t^O}  should  be  defined  as  a 
stochastic  process.  Our  first  requirement  Is  that  almost  all  sample 
paths  of  Y^  be  continuous.  This  eliminates  processes  obeying  the 
original  assumption  (A2)  (white  noise  process)  but  still  leaves  quite 
a large  class  to  choose  from.  To  retain  most  of  the  flavor  of  the 
origins''  fA2),  let  us  visualize  for  a while  that  the  search  region  A 
has  been  partitioned  Into  small  cells  AA  and  that  the  time  Interval 
I0,T]  has  also  been  divided  Into  small  subintervals  At,  and  let  us 
replace  (A2)  by: 

(A2 ' ) If  at  time  t the  observer  Is  In  a cell  AA 

then  at  t + At  he  la  equally  likely  to  be  In 
any  of  the  cells  adjacent  to  AA. 

In  other  words,  we  now  assume  that  the  observer  performs  a symmetric 
Bernoulli  random  walk  on  the  partition  of  A.  Regarding  this  as  a discrete 
approximation  to  a time-continuous  motion  in  a continuum,  we  obtain  our 
fundamental  assumption: 

A random  motion  is  to  be  modeled  as  a symmetric  (l.e., 

zero  drift)  Wiener  process. 

It  may  still  be  objected  that  no  moving  craft  can  in  reality 
follow  a sample  path  of  a Wiener  process.  This  is  true  since  the  sample 
paths,  although  being  continuous  with  probability  one,  are  almost  surely 
nowhere  differentiable.  Nevertheless,  we  are  at  least  closer  to  reality 
than  we  were  with  discontinuous  paths. 

In  the  remainder  of  this  section,  we  derive  a new  formula  for 
the  probability  of  detection  using  our  modified  model.  As  mentioned  in 
the  introduction,  we  restrict  ourselves  here  to  the  case  where  the  search 
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takes  place  In  a one-dimensional  region.  Note  that  although  the  basic 
formula  (1)  was  derived  In  [1]  for  planar  regions,  Its  derivation,  the 
problem  with  assumption  (A2) , and  In  fact  the  entire  discussion  so  far, 
applies  to  any  number  of  dimensions. 

For  one-dlmenslonal  search,  the  region  A will  be  an  Interval 
[-b,al,  a > 0,  b > 0.  We  consider  first  a stationary  target  ■ Xq 
uniformly  distributed  over  [-b,a]  with  the  observer's  location 
being  a Wiener  process  with  drift  p > 0 and  variance  parameter  > 0. 
Without  loss  of  generality,  we  choose  Yq  ■ 0,  the  observer  at  the 
origin  at  the  beginning  of  search.  Since  the  observer  clearly  should 
not  leave  the  search  region  [-b,a],  we  consider  the  enpolnts  -b,a 
to  act  as  reflecting  barriers.  Let  be  the  time  when  the  target 

1s  detected,  and  let  ■ min  {tSO;  Yj.“x}  be  the  first  time  the  process 
Y^  reaches  level  x.  With  c > 0 being  the  detection  distance  as  in 
(A3)  we  have  for  Xq  ■=  x the  relation 

T If  c < X s a, 

x-c 

0 If  |x|  ^ c,  (2) 

T ^ if  -b  i X < -c. 
x+c 


Typically,  we  assume  that  c <<  a+b.  Our  goal  Is  to  compute  the  proba- 
bility of  detection  by  the  time  t,  p=P(T^^^^t) , which  corresponds  to 
the  right-hand  side  of  (1)  with  L = vt.  Now  according  to  (2) 


-c 


P = 


(t) 


det 


1 

a+b 


P(T^_^«)dx  + ^ 


-b 


P(T^^St)dx, 


(3) 
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Let 


y(u,x) 


e"“‘^P(T  € dt) 

X 


be  Laplace  transform  of  the  distribution  of  T^.  Ry  taking  Laplace 
transform  of  the  diffusion  equation  for  the  Wiener  process  and  solving 
the  resulting  second-order  differential  equation,  we  obtain  (see  [2], 
p.  233  for  details) 


y(u,x) 


cosh  — 


cosh 

a 


cosh  — 
a 

cosh  ^ r/lu 
a 


if  0 < X i a. 


if  -b  i X < 0. 


(4) 


Hence  from  (3)  and  (4) 


a-c 


det 


Y(u,x)dx  + — 


Y(u,x)dx 


-b+c 


1 

a+b 


cosh  — /2i7 

a 


a-c 

r 


+ cosh  — /2u 

cosh  ° 


(5) 


b-c 


dx 


cosh  ^2u 
a 


Next,  using  the  expansion 


- — = 2 I z > 0 , 

n=0 


cosh  z 


8 


and 


cosh  — /in  - 4 
0 2 


+ e 


we  have 


a-c 


cosh  — 
a 


dx 


cosh  »'^2u 
0 


/lu 


i=lL 

2n+l 


n“0 


-(2n+l)  - /lu 


-2n  — /lu  -((2n+l)  (a+b-c)-b)  — /lu 


(6) 


+ e 


- e 


-((2n+l)  (a+b-c)+b)  /in 


- e 


and  similarly  for  the  second  term  In  (5).  However,  the  inverse  Laplace 


transform  to 

-K  /2u  , 

, K ^ 0,  is 

/hi  /H£ 


and  hence  inverting  (6)  we  obtain  for  the  density 
of  the  detection  time  the  expression: 


_d_ 

dt 


(t) 


(t) 

det 


2o 


1-  I 


(-1)° 


9 2 2 
2n  a 

to2 


(2n-l)(2n+l) 

n=i 


+ e 


2nV 

ta2 


-(2n+l)2(a+b-c)2  ^ 

2n+l  ® 


a2 

e'  cosh  (2n+l)aCa+b-cJ 

to2 


(7) 


+ e 
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Notice  that  the  infinite  aeries  inside  curly  brackets  converge  unifornly 

for  t 0,  and  very  rapidly.  Thus,  the  value  of  the  density  f„  (t) 

^det 

can  be  computed  approximately  for  each  fixed  t by  taking  only  the 
first  few  terms  of  the  series.  Rather  than  doing  that,  we  prefer  to 
look  at  the  asymptotic  behavior  of  the  density  for  small  t,  that  is 
at  the  beginning  of  search.  It  is  easily  seen  that  as  t (H* 


am. 


f».  (t) 

det 


F™  (t) 
det 


v8) 


(9) 


Thus,  at  the  beginning  of  search  the  probability  of  detection  by 
the  time  t Increases  like  vT,  in  contrast  to  the  corresponding 
asymptotic  behavior  obtained  from  (1),  where  the  Increase  is  linear. 

Remark!  For  a uniformly  distributed  stationary  target  the  distribution 

function  F,  ^(t)  can  also  be  derived  from  the  range  R^  of  the  Wiener 
det  t 

process.  The  range 

R^  = max  - min 
S£t  s^t 


is  the  length  of  the  interval  visited  by  the  process  Y^  up  to  the  time 
t.  Clearly 


for  c <<  a+b 


(10) 


10 


For  t >-0+  we  expect  E(R^)  to  be  almost  as  if  the  reflecting  barriers 
were  absent,  in  which  case  it  is  shown  in  [3]  (for  o ■ 1)  that 

E(R^)  “ so  that 

t 7T 

F.  ^(t)  ~ in  agreement  with  (9). 
det  a+b  Tt 

Next,  let  us  investigate  the  behavior  of  large  t,  that  is 

when  the  region  I-b,a]  is  becoming  saturated  with  search.  Performing 
the  integration  in  equation  (5)  we  have 

oo 

e (dt)  ” ^ cosh  — /lu 

^det  ^ 

(a+b-c)  /lu  b /2u 

arctan  + cosh  — /2u  (H) 

1 . (a+2b-c)  'flM  ^ 

1 + e 

(a+b-c)  a /2u 

e - e 

arctan  . 

, . (2a+b-c)  /2u 

Calling,  temporarily,  the  right-hand  side  of  (11)  w(u)  it  can  be  easily 
verified  that  for  each  u > 0 

uj(tu) 

; -►I  as  T -*•  0+  . 

cj(t) 

Hence,  by  Tauberlan  theorem  of  Feller  ([4],  p.  4A3) 

F (t)  ~ o)(^)  as  t 
det 


->•  CO 
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Since,  for  instance  w(^)  ~ e 


-(a+b) 


/! 


we  can  write  F_,  (t)  c 

^det 


-(a+b) 


fi 


as  t ' “ , 


which  is,  however,  not  very  informative  since,  in  fact. 


► a 


w(-)  ~ e 


1 - e 


- K 


for  any  a > 0 and  any  K. 

Let  us  now  compute  the  expected  value  of  the  detection  time.  From 


(2) 


ECT  ) 
x-c 


dx 

a+b 


a-c 

» 

. 


E(T  ) %-  +• 
X a+-b 


-b+-c 


E(T  ) 

X 


dx 

a+b  • 


(12) 


Since 


E(T  ) 

X 


- 3ir 


we  obtain 


— (x+2b)  if  0 < X i a , 


E(T  ) 

X 


— 2 (x-2a)  if  -b  i X < 0 . 


(13) 


Substituting  into  (12)  and  integrating  gives  the  result 


12 


- 7T 


1 


3a  (a+b) 
3 


(a-c)^  + 3b(a-c)^  + 3a(b-c)^ 


(lA) 


+ (b-c) 


* for  c <<  a+b  . 

3a 


Thus,  the  mean  detection  time  increases  approximately  as  a square  of  the 
length  of  search  region.  Higher  moments  of  the  detection  time  can  be 
obtain  d in  similar  fashion. 


13 


3.  Moving  Target 

In  the  previous  section,  the  target  was  assumed  to  be  stationary. 
Here,  we  would  like  to  consider  the  case  when  the  target  Is  also  under- 
going random  motion  over  the  search  region.  An  obvious  extension  of  the 
model  discussed  previously  would  be  Just  to  add  the  assumption  that  the 

target's  position  X^,  t ^ 0 Is  also  a Wiener  process  with  drift 

2 

“ 0 and  variance  parameter  > 0.  The  Initial  position  Xq  could 

still  be  assumed  uniformly  distributed  over  the  search  region. 

Since  by  (A3')  detection  occurs  as  soon  as  ^ ^ since 

■ X^  - 1s  again  a Wiener  process  with  drift  = 0,  variance 
2 2 2 

parameter  " °X  °Y  ’ Initial  distribution  same  as  that  of  Xq 

(we  still  assume  Yq  ■ 0) , the  detection  time  is  now  simply  the 

time  the  process  first  enters  the  interval  {-c,c].  Unfortunately, 

the  presence  of  a reflecting  barrier  at  -b  and  a in  the  original  model 

results  in  rather  complicated  boundary  conditions. 

We  will,  therefore,  restrict  ourselves  to  the  case  when  there  are 

no  reflecting  barriers.  For  Instance,  we  may  assume  that  the  initial 

distribution  of  Xq  Is  uniform  over  some  interval  [-bQ,aQ]  such  that 

3q  <<  a and  bQ  <<  b.  For  small  t the  distribution  of  detection  time 

should  be  approximately  the  same  as  If  the  barriers  were  absent. 

Let  be  the  first  time  a Wiener  process  W^  with,  Wq  = 0, 

2 

drift  p * 0 and  variance  parameter  o >0  crosses  level  x.  As  is 
well  known,  the  density  of  is  inverse  Gaussian,  namely, 

x2 

f_  (t)  = e ^°*',t>0,  - “><x<+“.  (15) 
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Since  for  Zq  ■ Zq  we  have  again 


det 


0 if  IzqI  ^ c , 


^Zq+c  ■*’0  ^ ^0  * 


we  can  write  for  the  density  of  the  detection  time 


^det 


Sq-C 


l/zTTO^t' 


2o^t  . 
e dx 


2 

2o  t . 
e dx 


2 3 

V2TTO  t 


2a 


l.i[r  2»'t  2a^t  I 


2 


2o„ 


where  o * a„ 


^0^+0^  . Notice  that  for  t -*■  0+  F_  (t)  ~ — 7^  )/— 
X Y Tdet  aQ'^^’o  ^ 


as  for  a stationary  target,  the  only  difference  being  in  the  variance 
parameter . 

Next,  consider  the  case  when  the  search  region  is  the  entire  line 

(no  reflecting  barriers)  and  the  Initial  distribution  of  Xq  is  normal 
2 

W(0,Ta^),  T > 0.  This  model  may  correspond  to  the  case  where  the  searcher 
arrives  at  time  t = 0 at  the  place  where  the  target  was  x time  units 
before.  The  density  is  now 
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However,  in  the  absence  of  reflecting  barriers,  tlie  mean  detection  time 

E(Tdet)  is  infinite,  in  fact,  for  any  initial  distribution  of  Xq  since 

the  inverse  Gaussian  distribution  (15)  has  infinite  mean  for  all  x 0. 

In  view  of  this  fact,  it  may  be  interesting  to  ask  whether  there  is 

continuous  motion  of  the  searcher  with  bounded  speed  which  would  result 

in  a finite  mean  detection  time.  More  precisely,  we  assume  that  the 

target's  motion  is  still  a Wiener  process  X^,  t ^ 0 with  drift  u = 0, 

2 

variance  parameter  ^ initial  distribution  Xq  normal 

2 

W(0  The  search  region  is  the  entire  real  line,  and  the  detection 

mode  as  in  (A3').  However,  we  now  allow  the  observer  to  choose  any  motion 
, t 2:  0,  deterministic  or  random,  with  = 0 as  long  as  his  speed 
between  turns  remains  bounded,  i.e.,  ^ v whenever  the  derivative 

exists. 


amn 
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It  seems  reasonable  to  assume  that  the  searcher  may  wish  to  move 
with  maximum  speed  v between  turns.  At  the  same  time  he  must  plan  on 
turning  Indefinitely  since  If  he  ever  (before  detection)  stopped  turning 
and  proceeded  straight,  the  mean  detection  time  would  still  be  Infinite. 
Hence,  we  will  assume  that  he  performs  a zigzag  motion  sweeping  with 
speed  V around  0 with  Increasing  amplitude.  A typical  path  of 

this  kind  of  motion  Is  depicted  in  Figure  1. 


FIGURE  1. 


Here  tj^,  t2>  . . . , 
Clearly,  with  tQ  ■■ 


are  turning  times  and  a = |Y  | 

n tn 

‘o  ■ “n  ■ 'n+1  - 


are  sweep  amplitudes. 
. . . , we  have 


vAt 

n 


n+1 


+ a 


(16) 


We  now  show,  following  a suggestion  of  A.  Washburn,  that  with 


n 

a = a , n 
n 


1,  2, 


> 


a > 1 , 


(17) 
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we  have  Indeed  *■  begin,  let  us  write 


■=<'^dctiv-)  -7 

o/t  \a/T  / 


+ 1 ^<'’detlV">  ■ 

o/t  \o/t  / 


whe  re 


(p(x)  = e ^ 


‘=C'detlV’‘>  ^ *(7=) 

o/t  \o/t  / 


-7  ♦(7=) 

o/t  \o/t  / 


M ^ ^ ^("det^‘^2n-ll 
0 L 


t , ” 


U '■0‘t,  7“2„-llVf>“2n-l  -7  ♦(7=) 

n=l  2n-l  o/t  \o/t  / 


P(X  >a„  , X_=x) . Next,  for  a ^ x 

4n-l  2n-l'  0 


P(X|.>a|XQ=x)  = 


J 0/7  ‘ 

\o/7, 

1 

!)/— 

1 y” 
o/t 

Ao/T 

d)/ 

^a-x  \ 

a-x  1 

ko/t"/ 
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so  that 


dx 


2 ^2n-l 


■Szl_  7!2S=iZ\!!2l!zi  JJL.) 
'2n-r*  /r  Vo/^/ 


dx 


2 ^ 2n-l 


o/7  Va/r/  \0'^t2n_i+T  J 


2 ^ 2n-l 


At 


2n-l 


/? 


2n-l 


^2n-r^ 


X - 


2n 


-i\ 


T+t 


2n-l 


\>'-2„-l  / 


dx  + At 


2n-l 


a/7 

(0/7) 

dx. 


2 ® 2n-l 


Now 


2 ® 2n-l 


-3-  J-3S-\  dx  X.  ♦(— ) 

a/T  \a/x  J ^2n-l  \2o/r  / 
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and  since  0<x<-;ra„  ,■» < < 


2 2n-l  a2|^_j_  ^2n-l  ^ ®2n-l 


2 ® 2n-l 


/ 


/E 


2n-l 


^2n-l"’' 


x - 


2n-l 


T+t 


2n-l 


J *^2n-l 

\ / 


2n-l  -1/2- 


2 raax{l,T  } . 


2n-l 


Substituting  into  (19)  we  obtain  the  bound 
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) 


ao 


I ^«detlV> 
0 


1 

o/r 


Cl 

dx  S ^ + I 


1 


,,  -1/2,  ^*^2n-l  *^-1 

2 raaxll.T  ) 

^2n-l 


. I 

n-l  “2n-l  \,2o/7  / 

Now  with  as  In  (17),  the  latt<;r  series  obviously  converges  and  so 

does  the  former,  since  by  (16) 


^•^2n-l  ‘^‘^2n-l  n.  ^2n-l  n, 

= 0(a  ) , ■ ■ ■ * 0(a  ) , 

2n-J 

” 2n 

and  'l  a e is  a convergent  series.  Hence,  the  first  Integral 

n=l 

In  (18)  Is  finite.  The  second  integral  is  handled  in  the  same  fashion. 

It  would  be  of  considerable  interest  to  determine  the  motion  of 
the  searcher,  constrained  to  be  continuous  and  of  bounded  speed,  which 
actually  minimizes  the  mean  detection  time.  Although  one  may  conjecture 
that  the  optimum  path  will  be  as  in  Figure  1,  the  optimal  choice  of  the 
turning  times  t^  remains  an  open  problem. 
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